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ABSTRACT 

We present B, V, I time-series photometry of the Coma Berenices dwarf 
spheroidal galaxy, a faint Milky Way satellite, recently discovered by the Sloan 
Digital Sky Survey. We have obtained V,B — V and V,V — I color-magnitude 
diagrams that reach V ~ 23.0 — 23.2 mag showing the galaxy turnoff at ~ 21.7 
mag, and have performed the first study of the variable star population of this 
new Milky Way companion. Two RR Lyrae stars (a fundamental-mode -RRab- 
and a first overtone -RRc- pulsator) and a short period variable with period 
P=0. 12468 days were identified in the galaxy. The RRab star has a rather 
long period of Pat = 0.66971 days and is about 0.2 mag brighter than the RRc 
variable and other non-variable stars on the galaxy horizontal branch. In the 
period-amplitude diagram the RRab variable falls closer to the loci of Oosterhoff 
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type-II systems and evolved fundamental-mode RR Lyrae stars in the Galactic 
globular cluster M3. The average apparent magnitude of the galaxy horizontal 
branch, (Vhb) = 18.64±0.04 mag, leads to a distance modulus for the Coma dSph 
/io = 18.13 ± 0.08 mag, corresponding to a distance d=42'^ kpc, by adopting a 
reddening E{B -V) = 0.045 ±0.015 mag and a metallicity [Fe/H] =-2.53 ± 0.05 
dex. 

Subject headings: galaxies: dwarf — galaxies: individual (Coma) — stars: dis- 
tances — stars: variables: other — techniques: photometric 



Introduction 



Over the past few years the analysis of the Sloan Digital Sky Survey (SDSS) data 
led to the discovery of several ultra-faint companions of the Milky Way ( MW). The new 
syste ms include two very faint globular clusters (GCs): Koposov 1 and 2 (IKoposov et al. 



20071 ). and fourteen dwarf spheroidal (dSph) galaxies: WillmanI, Ursa Major I, Canes Ve- 



Venaticill fCVnII), Leo IV, 


iercules. 


1/eo T, Bootes II, Leo V, and Bootes III ( Wi 


— ^, ^^^^^^ 
Iman et al. 
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t al. 2006a: 


Grillmair 


2006: Zucker et al.ll2006b: Belo 


curov et al. 


2006. 2007: 


Irwin et al. 


2007; 


Walsh et al. 2007: I 


3elokurov et al.l2008: Grillmaii 


2008). The 


new dSphs 



have half-light radii similar to those of the classical MW dSph companions, however, they 
are much fainter, with typical effective surface brightnesses in the range from 28 to 31 



arcsec 



( iBelokurov et al.ll2007l . hereafter B07). The new discoveries along with the 10 pre- 
viously known MW dSph companions (Draco, Ursa Minor, Fornax, Carina, Sculptor, Leo I, 
Leo II, Sextans, Sa gittarius and Canis Major; lMatedll998l : Ibata, Gilmore, & Irwin 1995; 
Martin et al.ll2004l ). bring to twenty- four the number of dSph galaxies presently known to 
surround the MW. In the absolute magnitude versus half-light radius plane the new dSphs 
are well separated from both GCs and traditional dSph galaxies, since they are generally 
more extended than GCs and much fainter than the old dSphs (see Fig. 8 of B07). Only 
CVnl, the brightest of the new SDSS dSphs, lies on the faintest tail of the distribution for 
traditional dSphs. All the new SDSS dSphs host a metal-poor and ancient stellar compo- 
nent j_with_jii^taJ]ici^^ —3 
dex JMuuoz et al. Ibood : Isimon fc Gehall2007l : ICeha et al Eooi . iKirbv et al. ]l2008h . and age 
as old as the stars in Galactic GCs like M92 (see B07). The density contours are irregular 
since most of these galaxies are in process of disruption. The closest ones ( Bootes I, UMall 



and Coma) are more distorted o wing to the tidal interaction with the MW (IBelokurov et al. 



2006. 


2007: 


Zucker et al. 


2006b) 
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the exception only of Leo T, the new systems all have globular-cluster-like CMDs, showing 
main sequences, turn-offs, as well as hints of red giant and horizontal branches. Variable 
stars have been i dentified so far i n thr e e of the new systems, nam ely Bootes I, CVn I, and 
CVnII. Bootes! JPaU'Ora et all bood : Isiegell bood ) and CVn II JCreco et ali boosh both 
were found to con tain RR Lyrae stars with properties resembling those of Oosterhoff-type II 



Oosterhofl 



19391) G Cs in the MW. Confirming its similarity to the traditional dSphs, CVn I 



( iKuehn et al.ll2008l ) hosts, instead, RR Lyrae stars with properties intermediate between the 



two Oosterhoff types, along with a few candidate Anomalous Cepheids. 

The Coma galaxy {R.A. = 12^'26°>59.0^ DEC = +23°54'15", J2000.0, Z=241.9°, 6=83.6° 
B07), located at a heliocentric distance of 44 ± 4 kpc with a position angle of 120° and an 
absolute magnitude of My = —3.6 ±0.6 mag, is among the faintest of the SDSS new discov- 
eries. It has half light radii: r ^{Plummer) = 5.0 arcmin, and rh{exponential) = 5.9 arcmin, 
respectively, and an irregular and extended shape (B07). The i, g — i CMD obtained by B07 
goes down to i ~ 25 mag, is consistent with a single, old stellar population of metallicity 
[Fe/H] ~ —2 , and shows a poorly populated horizon tal branch (HB) . Meta llicities of radial 
velocity members of the Coma dSph were obtained by lSimon fc Gehal (120071 ) an d confirm the 



l ow m etal abundance of the galaxy. These metallicities were re-evaluated by iKirby et al. 



( 120081 ) ■ who derived for the galaxy an average metallicity ([Fe/H]) = —2.53 ± 0.05 dex with 
dispersion o"[Fe/H]=0.45 dex. 

In this Letter we present the first study of the variable stars in the Coma dSph, which led 
us to identify two RR Lyrae stars and one short period variable in the galaxy. We provide 
multiband light curves and properties of the variables, and use the magnitude of the HB 
stars to measure the distance to Coma. We also present B, B — V and V, V — I CMDs to 
V ~ 23.0 — 23.2 mag, which show the galaxy main sequence turnoff at ~ 21.7 mag, and 
reveal the presence of several candidate blue straggler stars (BSSs). 



2. Observations and Data Reduction 



Time-series B,V,I photometry of the Coma dSph galaxy was collected on 2007, March 
12-16 and April 13-17, with BFOSC at the 1.52 m telescope of the Bologna Observatory 
in Loianc0 {V and / data), on April 30 a nd May 9-10, using WIRO-Prime, the prime 
focus CCD camera (jPierce fc Nations! |2002| ) of the 2.3 m Wyoming Infrared Observatory 
telescope (WIRO) {B and V data), and B,V,I observations were gathered on April 21- 
22, using WFC (Wide Field Camera), the prime focus mosaic CCD camera of the 2.5 m 
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Isaac Newton Telescope (INT). The field of view (FOV) covered by the three instruments 
is 13 X 12.6 arcmin^ for BFOSCOLOIANO, 17.8 x 17.8 arcmin^ for WIRO, and 33 x 33 
arcmin^ for WFC@INT. We needed two partially overlapping fields to cover the galaxy with 
BFOSC, while just one WIRO pointing was sufficient, and from the WFC@INT field we 
could also infer some additional information about the contamination by field stars and 
background galaxies. We obtained 14 B images, 83 V frames and 64 / frames of the Coma 
dSph. Images were pre-reduced following standard procedures (bias subtraction and fiat- 
field correction) with IRAF. The /-band images were corrected for fringing by adopting, for 
each instrument, a well-suited fringing map. We the n perfornaed PSF ph otometry using the 
DAOPHOTIV/ALLSTAR/ALLFRAME packages JstetsonI EqStI . Il994h . Typical internal 
errors of the single-frame photometry for stars at the HE magnitude level are of about 
0.005 mag in all the three bands. The absolute photometric calibration was derived using 
observations o f standard stars in the Landolt fields SA 101, SA 107, SA 110 and PG1323 
jLandoltl[l992l ). as extended by P.B. Stetsoifl, which were obtained at the INT during the 
night of 2007, April 22. Errors of the absolute photometric calibration are as = 0.01, 
av = 0.01 and aj = 0.02 mag, respectively. 



3. Identification of the Variable Stars 



Variable stars were identified using t he V and / time-series data, separately. First we 
calculated the Fourier transforms (in the I S chwarzenber g- Czernyl Il996l formulation) of the 
stars having at least 12 measurements in each photometric band, then we averaged these 
transforms to estimate the noise and calculated the signal-to-noise ratios. Results from the 
V and I photometries were cross-correlated, and all stars with S /N> 4 in both photometric 
bands were visually inspected, for a total of about 500 candidates. We also checked whether 
some of the stars in the BSSs region might be pulsating variables of SX Phoenicis (SX 
Phe) type. Study of the light curves and period derivation were carried out using GRaTiS 
(Graphical Analyzer of Time Series, IClementini et al.ll2000l ). We confirmed the variability 
and obtained reliable periods and light curves for 2 RR Lyrae stars: 1 fundamental-mode 
(RRab) variable with period P = 0.66971 days (VI), 1 first-overtone (RRc) pulsator with 
period P = 0.31964 days (V2), and for one short period variable with period P = 0.12468 
days (V3). Identification and properties of the confirmed variable stars are summarized in 
Table [H their light curves are shown in Fig [TJ Time-series data for the variable stars are 
available on request from the first author. Both the short period variable and the RRc star 
lie inside the half-light radius of the Coma galaxy. The RRab lies, instead, slightly outside 
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this region. In the CMD, the short period variable is located in the region of the BSSs, thus 
suggesting it might be a n SX Phe star. However, the star does not fit t he SX Phe period- 
luminosity relations (e.g. ISantolamazza et al. II2OOII . iPoretti et al. II2008I ) unless a pulsation 



period of one third the value supported by the period search analysis is adoped. The shape 
of the light curve and the similarity of the Ay, Aj amplitudes may suggest a classification as 
a W UMa binary system. However, the present data do not allow us to definitely assess the 
actual nature of the star. 

The two RR Lyrae stars both fall near the HB of the Coma dSph. However, while 
the average magnitude of the RRc star is fully consistent with the average luminosity of 
the HB we infer, by fitting the galaxy CMD with the ridge lines of the Galactic GC M 68 
(see Section 4), that the RRab star is about 0.2 mag brighter. The overluminosity of the 
RRab star could be due to geometrical projection effe cts, since the galaxy appears to be 
elongated due to the tidal interaction with the MW (see lDalPOra et al.ll2008l and Section 4). 
Contamination by an unresolved companion may also be possible, but this is constrained by 
the large amplitude of the RRab star, and by its color, which is consistent with the expected 
color of an RRab star. A further possibility is that evolution off the zero-age HB contributes 
to the enhancement of the brightness of the RRab star (see end of Section 4). 



The rather long period of the RRab star suggests that Coma is, like CVn II (iGreco et al. 



20081 ). more similar to Oosterhoff type-II systems in the MW, and diff ers from CVn I, th e 
brightest of the SDSS dSph's, which has Oosterhoff- intermediate type (IKuehn et al.ll2008l ). 
Fig [2] shows the position of the Coma RR Lyrae stars in the l^-band period-amplitude 
diagram of the Bootes I and CVn II RR L yrae stars. The Coma RRab star (filled circle) 
falls closer to the loci of OoII systems (from lClement fc Rowell2000l . dashed- dotted line) and 



evolv ed fundamental-mode RR Lyrae stars in the Galactic GC M3 (from ICacciari et al. 



2OO5I ). The RRc star falls at the short extreme of the bell-shaped distribution defined by 
OoII RRc stars. Both these results suggest an OoII classification for the Coma dSph, 
although with some caution given the very small number of variable stars. We have used 
the parameters of the Fourier decomposition of the l^-band light curves, along with Jurcsik 
& Kovacs (1996) and [Morgan, et al. I (120071) formulas, to e stimate individual metallicities 
for the Coma RR Lyrae stars, on the IZinn &: WestI (119841 ) metallicity scale. Results are 
summarized in the last column of Table 1, they confirm the low metallicity derived for the 
Coma dSph by the spectroscopic studies. 
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4. CMD, Structure and distance 



The upper panels of Fig. [3]sliow the V, B — V and V, V — I CMDs of the Coma dSph. 
We have plotted all stellar-like objects (selected using the CHI and SHARP parameters 
provided by the reduction packages), over the whole 0.5 square degrees covered by the INT 
observations. This selection is rather secure for magnitudes above V ~ 20.5-21.0 mag, while 
it becomes increasingly uncertain at fainter magnitudes. The variable stars (filled triangles) 
are plotted according to their intensity-averaged magnitudes and colors (see Table 1). The 
CMDs reach V ~ 23.2 mag, and appear to be heavily contaminated at every magnitude level 
by field objects belonging to the MW halo and disk, and by galaxies which escaped selection. 
We have used the mean ridge lines of the Galactic globular cluster M68 (NGC4590; solid 
blue lines in Fig. |3]), drawn from Walker's (1994) B, V, I photometry for the cluster, shifted 
by AV^ = +2.95 mag in magnitude, and A{B -V) = +0.025 mag, A{V - I) = +0.030 mag 
in color, to match the horizontal and red giant branches of the Coma dSph. Then we selected 
as most likely belonging to the galaxy the sources lying within ±0.06 mag in B — V, and 
±0.05 mag inV — I, from the ridge lines of M68 (red dots in the upper panels of Fig. [3]). The 
slightly larger range adopted for the B — V color takes into account the larger uncertainty 
caused by the small number of B frames. To allow for the larger photometric errors below 
V = 21.5 mag, we also considered as belonging to the galaxy stars with colors in the range 
from ± 0.05/0.06 to ± 0.1 mag from the ridge lines of M68 (cyan dots). M68 is well suited 
for identif ying members of the Coma dSph since, like Coma, is very metal poor ([Fe/H]M68= 
—2.1 dex, I Walker Ill994l ) and has a well defined and tight red giant branch, as w ell as an 



exten ded HB well populated inside and to the blue of the RR Lyrae instabilitv strip (jWalker 
1994h . Adopting for M68 the reddening value of E{B - V)=0.07 ± 0.01 mag JWalker Ill994 ) 



the color shifts needed to match the horizontal and red giant branches of Coma imply a 
reddening of E{B — V) = 0.045 ±0.015 mag f or the galaxy, which c an be compared with the 
value of 0.019 ± 0.026 mag derived from the ISchlegel et al.l (119981 ) ma ps. We also show in 
Fig. 3 results of the fit b ased on a less metal poor globular cluster: M3 (IFerraro et al. 1119971 . 
Johnson fc Bolte Ill998l ). The red giant branch of M 3 is too red and, to match that of Coma, 
would require for the galaxy a negative reddening E{B — V) ^ —0.06 mag. With the help 
of the M 68 ridgelines it is possible to determine the average luminosity of the Coma dSph 
HB roughly at the center of the RR Lyrae instability strip: (liiB) =18.64 ± 0.04 mag, and 
to locate the galaxy main sequence turn-off at ^ ~ 21.7 mag. The CMDs also show several 
stars in the BSS region (green dots in Fig. [3]), which are likely BSSs of the Coma dSph. Our 
identific ation of the C o ma ra ember stars is confirmed in an excellent way by the comparison 
with the iKirby et al. I (120081 ) spectroscopic study of 24 stars with membership to the Coma 
dSph confirmed by radial velocity measurements (open circles in Fig. [3]). The agreement 
between the two studies is impressive, and proves the reliability of our procedure to select 



- 7- 



member stars, and our identification of tlie Coma HB. Tlie lower panels of Fig. [3] show the 
position of the stars we consider to belong to the Coma galaxy in the FOV of the 5-band (left 
panel) and \^-band (right panel) INT observations. Symbols and color-coding are the same 
as in the top-panels of the figure. Symbol-sizes are proportional to the object's magnitudes. 
The HB stars are outlined by squares. It is noteworthy that, with the exception only of the 
RRc variable, all the HB stars fall outside the galaxy half-light radius, but generally close 
to stars whose membership to Coma is confirmed by the spectroscopic studies. 

A further feature is visible in the CMDs of the Coma dSph (see upper panels of Fig. [3]), 
formed by sources having roughly 20.6 < V <21.9 mag, 1.3 < V — I < 1.5 mag, and 
B — V ~ 1.5 mag. Among them those marked by black dots are mainly concentrated in 
the left-hand portion of the maps outhned by dashed circles (see lower panels of Fig. [3]). In 
the literature, at this position, we find the cluster of galaxies MaxBCG J186. 85861+2380004 
( iKoester et al.l 120071 ). which has a photometric redshift z = 0.243050. The feature we see 
in the CMDs is, probably, the "red sequence" of this cluster, as confirmed by the V — I, 
B — V colors we derive for early-type galaxies at z~0.25 using the GISSEL03 spectral code 
jBruzual fc Chariot l^W^ . 

The average magnitude of the Coma dSph HB, set by the matching with the ridge lines of 
M 68, is (^hb) = 18.64±0.04 mag (where the error includes both uncertainties in the photom- 
etry and in the matching procedure). Assuming Mv=0.59=fc0.03 ma g for the absolute magni- 
tude of HB stars at [Fe/H]=-1.5 dex JCacciari fc Clementini|[2003h . AMv/A[Fe/H]=0.2U 
(=b 0.047) mag/dex for the slope of the luminosity-met allicity relation of the RR Lyrae stars 



(I Clement ini et al. 



dex ( Kirby et al. 



2003h . E{B - 1/)coma=0.045 ±0.015 mag, and [Fe/H]coma=-2.53 ± 0.05 



2008i) . the distance modulus of Coma is 18.13 ± 0.08 mag, ((i=42l^ kpc). 



Errors include uncertainties in the photometry, reddening, metallicity, and RR Lyrae abso- 
lute magnitude. On the same assumptions, the average apparent magnitude of the RRab 
star {V) = 18.44 ± 0.03 mag leads, instead, to the distance modulus of 17.93 ± 0.08 mag, 
{d=39'tl kpc). While the first estimate is in good agreement with the distance of 44 ± 4 kpc 
found by B07, the value inferred from the RRab stars is, respectively, 3 and 5 kpc shorter 
than the galaxy distance derived in the present study and in B07. An additional distance 
estimate for the two RR Lyrae stars was obtained using the empirical Wesenheit relation in 
the B, V bands as defined by Kovacs and Walker (2001) and calibrated by Kovacs (2003, 
this calibration is consistent with a distance modulus for the Large Magellanic Cloud of 
18.55 mag). The resulting moduh: /io(RRab) = 18.06±0.10 and /io(RRc) = 18.09±0.10 mag, 
are in good agreement with each other and are consistent with the distance inferred from 
the average magnitude of the galaxy HB. The overluminosity of the RRab reflects into the 
Wesenheit index, however it is balanced by the star longer period, thus the resulting dis- 
tance modulus is in good agreement with that derived for the RRc variable. These results 
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suggest that the RRab star is brighter because it has evolved off the zero age HB, rather 
than because it is significantly closer to us. 



5. Summary and conclusions 

We have identified and obtained V, I light curves for 2 RR Lyrae (1 RRc and 1 RRab) 
and one short period variable, in the Coma dSph galaxy. The behavior of the RR Lyrae stars 
suggests an Oosterhoff II classification for the galaxy. From the average luminosity of the 
HB stars, the distance modulus of the Coma dSph is /io=18.13 ± 0.08 mag {d = 42t^ kpc). 
The galaxy appears to be elongated with irregular and extended shape, likely caused by the 
tidal interaction with the MW. The RRab star lies outside the galaxy half-light radius and 
is about 0.2 mag brighter than the average magnitude of the galaxy's HB. Its bright mean 
magnitude leads to a distance of d=39 kpc, which, if interpreted in terms of projection 
effects, could provide a rough estimate for the galaxy elongation. However, the Wesenheit 
results suggest that the star overluminosity is indeed due, at least in part, to evolution off 
the zero-age HB. Finally, from the analysis of the CMDs we identify in the lower north-east 
f ield of the Coma ga laxy, the background cluster of galaxies MaxBCG J186. 85861+2380004 



( IKoester et all 120071 ). 
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Table 1: Identification and ])ro])erties of varial)le stars in tlie Coma dS\)li galaxy 



Name 


a 


S 


Type 


P 


Epoch(max) 


<V) 




Av 




[Fe/H]zw 




(2000) 


(2000) 




(days) 


(-2450000) 


(mag) 


(mag) 


(mag) 


(mag) 


a 


VI 
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"^Metallicities derived from the Fourier parameters of the hght curves. 
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Fig. 1. — V (left panels) and I (right panels) light curves of variable stars discovered in 
the Coma dSph. From top to bottom: fundamental-mode RR Lyrae star, first-overtone RR 
Lyrae pulsator, short period variable. 
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Fig. 2. — Period- amplitude diagram in the V band. D ashed and dotted-dashed hnes are 
the positions of the Ool and OoII Galactic GCs, from I Clement fc Rowd (120001 ). Period- 
amplitude distributions of the bona fide regular ( solid curve) and well-ev olved (dashed curve) 
fundamental-mode RR Lyrae stars in M3, from ICacciari et al. I (120051 ). are also shown for 
comparison. 
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Fig. 3. — Upper panels: V, B — V (left) and V, V — I (right) CMDs of sources in the INT 
field of view. Solid (blue) and dashed (black) lines are the ridgelines of the Galactic GCs 
M68 and M3, respectively. Red and cyan dots are stars respectively within ±0.06 mag in 
B — V and ±0.05 mag in V — I, and from ±0.05/0.06 to ±0.1 mag from the ridgelines of M68. 
Blue dots are non- variable stars on the HB, green dots are stars in the BSSs region. Variable 
stars are marked by filled triangles, cyan: RRc star, magenta: RRab st ar, blue: short perio d 



variable. Open circles mark 24 member stars of the Coma dSph from iKirby et al. I (120081 ). 



Black dots are galaxies of cluster MaxBCG J186.85861±2380004. Lower panels: maps of 
objects most likely belonging to the Coma dSph, color-coded and with same symbols as in 
the top-panels of the figure. HB stars are outlined by squares. The large circles show the 
galaxy half-light radius rh{exponential)= 5.9 arcmin (B07). The dashed circles show the 
position of the galaxy cluster. 



